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We perform Monte Carlo investigations of the 4d jV = 1 supersymmetric Yang-Mills (SYM) the- 
ory on the lattice with dynamical gluinos in the adjoint representation of the SU(2) gauge group. 
Our aim is to determine the mass spectrum of the low-lying bound states which is expected to be 
organised in supermultiplets in the infinite volume continuum limit. For this purpose we perform 
simulations on large lattices, up to an extension L/ro ~ 6 where ro ~ 0.5 fm is the Sommer scale 
parameter. We apply improved lattice actions: tree-level improved Symanzik (tlSym) gauge ac- 
tion and in the later runs a Stout-smeared Wilson fermion action. The gauge configuration samples 
are prepared by the Two-Step Polynomial Hybrid Monte Carlo (TS-PHMC) update algorithm. 
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1. Introduction 

The jV = 1 supersymmetric Yang-Mills (SYM) theory is the minimal SUSY extension of the 
SU(/V C ) Yang-Mills theory. The fermionic degrees of freedom, the gluinos, are given by the super- 
partners of the gauge fields (gluons) and are described by Majorana spinors X a (a = 1 . . . N% — 1) 
transforming according to the adjoint representation of the gauge group. SYM is characterized by 
a rich low-energy dynamics with interesting aspects as confinement and spontaneous breaking of a 
discrete chiral symmetry - a continuum U(l) chiral symmetry is missing at the quantum level due 
to the Adler-Bell-Jackiw anomaly. SYM can be related to QCD with a single quark flavour (Nf = 1 
QCD), where the Majorana spinor is replaced by the single Dirac spinor. The latter model is also 
object of investigation by our collaboration [jl|]. 

This work represents a continuation of a long-standing project of the DESY-Munster-Roma 
Collaboration (DMRC) for the simulation of SU(2) SYM, see [Q] for a review and [J3j] for more 
recent results. Following [|[] we apply the Wilson approach, which has been proved to be successful 
in lattice QCD computations in spite of its known limitations. SUSY is broken by the lattice 
discretisation and, in the Wilson approach, by the Wilson term. It is expected to be recovered 
in the continuum limit by properly tuning the only relevant parameter, the bare gluino mass, to a 
critical value corresponding to massless gluinos. Another (related) inconvenience of the Wilson 
discretisation, namely a non positive-definite fermion measure even for positive gluino masses, 
turns out to have no appreciable impact in practical applications. 

Past simulations of DMRC were performed on quite fine lattices, 1 a ~ 0.08 fm, but in a small 
volume, L ~ lfm; this setup was appropriate for the study of the SUSY Ward identities also 
valid in a finite volume. We now concentrate on the mass spectrum of bound states for which low- 
energy effective theories predict a reorganisation of the masses in two supermultiplets at the SUSY 
point []6[ ^]. Thanks to a new more efficient simulation algorithm (see below) and enhanced com- 
puting resources we are now able to accumulate relevant statistics on larger lattices. Our present 
series of numerical simulations are performed on 16 3 • 32 and 24 3 -48 lattices with lattice spacing 
a ~ 0. 125 fm. The lattice extension L ~ 2 — 3fm is expected to be large enough to allow control 
over finite volume effects on the bound states masses. Simulations on finer lattices are planned. 

2. Lattice formulation and algorithms 

2.1 Lattice formulation 

We apply the tree-level Symanzik (tlSym) improved gauge action for the gauge part including 
rectangular Wilson loops of perimeter six: 

S g = j3(c £{l - -UeTrC/ p/ } + Cl £{1 - -^ReTrt/ rec }) , (2.1) 

pi ly c rec ly c 

with Co = 1 — 8ci and ci = — 1/12 in the case of tlSym action. 

The contribution of the gluino to the effective gauge action is given by 

Sf = -ilogdetg[t/] , (2.2) 
'We use QCD units for setting the scale. The Sommer scale parameter is fixed to the value ro = 0.5 fm. 
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where Q is the non-hermitian Dirac -Wilson fermion matrix defined by 

Q*[U] = 8^ b - K £ (4, y+ A (1 + rv)Vf(y) + 8 x+fl>y (l - Y»)V™(x)); (2.3) 
V^(jc) is the adjoint gauge link, a real orthogonal matrix: 

v£>[u] (,) = 2Tr{c/t W r^^r 6 } = y™ fo (x) = [y- lafc (*)] r . (2.4) 

(r a are the generators of the SU(/V C ) group. In case of SU(2) one has T a = \o a with the Pauli 
matrices a a .) Observe that det<2 > for Majorana fermions. In our recent simulations the links 



U x> ^ in the Wilson-Dirac operator Eq. ( |2.3[ ) are replaced by stout-smeared links [g], which are 
defined as 

U§ = U Xrtl exp | i (fl^ - fttj - i Tr (ft^ - QtJ | , ^ = U^C^ (2.5) 
where C v >fl is a sum over staples 

CjC,/i = V, PuV (Uy_i_h yUx+V.U.Ux,V -\-Ux-V+U.vUx-V.U.Ul-?i y) • (2-6) 

We apply one step of smearing with smearing parameters p^ v = 0.15 in all lattice directions. 

2.2 Updating algorithm 

The factor in front of logdet(<2) in the Curci-Veneziano action reproduces the absolute value 
of the Pfaffian for Majorana fermions. Effectively, it corresponds to a flavour number Nf = \. 
The sign of the Pfaffian can be included in a reweighting procedure (see below). In our numerical 
simulations we use the two-step polynomial hybrid Monte Carlo (TS-PHMC) algorithm [§] which 
is based on a two-step polynomial approximation [[K]]. The fermion determinant is represented as 

detiQff = {det( Q?Q)} N f' 2 ~ — t-t— (2.7) 



with the condition on the polynomials P n 



lim P n Ax)P n Jx) =x- N fi\ xe [e,X], (2.8) 

where the interval covers the spectrum of Q^Q. The order of the first polynomial n\ is chosen as low 
as possible provided that the acceptance in the accept-reject step done with P„ 2 , after a sequence of 
PHMC trajectories prepared with P nj , is sufficiently high (in our cases 80-90%). 

In earlier simulations of the DESY-Munster collaboration the two-step multi boson (TSMB) 



algorithm [JlOp was used for the SYM theory. There the updating is performed by heatbath and 
overrelaxation sweeps for the pseudofermions and Metropolis and heatbath sweeps for gauge field. 
In case of the TS-PHMC algorithm, which is more effective in producing short autocorrelations 
among the gauge configurations, PHMC trajectories are created by applying the Sexton- Weingarten 
integration scheme with multiple time scales and, as usual, Metropolis accept-reject at the end of 
each trajectory. After the update sweeps/trajectories the second high-precision polynomial P n2 is 
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used in a stochastic noisy correction step with a global accept-reject condition. The polynomial 
P ni , which is an approximation for [x~ N f/ 2 P ni ]~ l , has to be highly precise so that the error of the 
approximation is negligible compared to the overall statistical error, provided that the spectrum of 
Q}Q is in the interval [e, A]. 

If some of the eigenvalues is outside [e,A] then the approximations are not precise enough 
and a reweighting procedure has to be applied by an appropriately chosen high order polynomial. 
This gives a reweighting factor C[U]. In this reweighting procedure the sign of the Pfaffian of the 
fermion matrix can also be included in the measurement as follows: 

(signPf[£/]C[£/]A[£/]) g 
{ > (signPf[C/]C[C/])„ • 1 • J 



3. Simulation details 

The algorithmic parameters of our TS-PHMC runs are summarised in Table ||. The runs are 
performed at a single lattice spacing corresponding to j8 = 1.6. The measured Sommer parameter 
for the lighter gluino masses is about ro /a ~ 4. By using the QCD scale as mentioned in the 
Introduction we obtain a ~ 0.125fm and a box-size L ~ 2fm on the 16 3 • 32 lattices and L ~ 3fm 
on the 24 3 -48 lattices. Comparison of results from ensembles Ci a m on the 16 3 • 32 lattice and C on 
the 24 3 • 48 lattice at the same gluino mass allows the study of finite size effects. 

An indication of the lightness of the gluino is given by the adjoint pion mass, M a - n , extracted 
from the connected part of correlator of the pseudoscalar gluino bilinear (see below). The lightest 
simulated adjoint pion mass (run D) corresponds to M a - K ~ 353(20) MeV in QCD units. 

Compared with the previously used TSMB algorithm, TS-PHMC displays a substantially im- 
proved update efficiency with shorter plaquette autocorrelation times x plaq . This is particularly true 
in the light gluino regime where the efficiency of TSMB undergoes a strong depletion. 

Some runs required the computation of the correction factor and the determination of the 
Pfaffian sign. We computed the correction factors for all runs in Table [I] with the exception of runs 
A, B, A s and B s . Most signs of the Pfaffian are positive; the run with largest number of negative 
Pfaffians is run D where we found 15 configurations out of 5160 with negative sign. In most cases 
the effect of the correction factor turned out to be negligible. This was not the case in run D where 
the effect of C[U] was important for the masses of adjoint meson bound states and the gluino- 
glueballs. All quoted statistical errors of the measured quantities were estimated by the T-method 

n 



4. Bound states 



The bound states masses are computed from the zero-momentum correlation function of the 
corresponding interpolating operator 0. 

i) Adjoint mesons. Low-energy theories [Q, ^] predict a Wess-Zumino supermultiplet con- 
taining colourless composite states of two gluinos. Such states have spin-parity quantum numbers 
0~ and + . In analogy to flavour singlet QCD we denote the former a-r\' and the latter a-fo. To 
project these states on the lattice we used the gluino bilinear operators G = ATA where Y = 75, 1 
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Table 1: TS-PHMC runs algorithmic parameters with tlSym at j3 = 1 .6. Runs labelled with s have been per- 
formed with Stout-links; 8l refers to the total trajectory length and A^c is the noisy correction acceptance. 



Run 


L 3 .T 


k 


#Traj. 


St 


A N c% 






e 




A 


"l 


«2 


A 


16 3 .32 


0.1800 


2500 


1.05 


95.6 


7.5 


4.25 


•10 


-3 


3.4 


50 


100 


B 


16 3 .32 


0.1900 


2700 


1.05 


96.4 


3.08 


9.5 


10 


-4 


3.8 


80 


300 


La 


Id .32 


0.2000 


1 ATI 

1973 


0.99 


82.9 


4.91 


5.0 


10 


-5 


4.0 


200 


700 


Cb 


16 3 .32 


0.2000 


8874 


0.99 


88.3 


27.6 


5.0 


10 


-5 


4.0 


200 


700 


D 


16 3 .32 


0.2020 


6947 


0.56 


88.5 


45.7 


1.0 


10" 


-6 


4.0 


800 


2700 


A 


24 3 .48 


0.1980 


1480 


0.9 


89.6 


4.64 


1.0 


10 


-4 


4.0 


200 


600 


B 


24 3 .48 


0.1990 


1400 


0.9 


88.7 


2.65 


4.0 


10 


-5 


4.0 


270 


800 


C 


24 3 .48 


0.2000 


6465 


1.0 


88.6 


7.4 


2.0 


10" 


-5 


4.0 


350 


1000 


A s 


24 3 .48 


0.1500 


370 


1.0 


97.3 


3.5 


5.5 


10 


-5 


2.2 


200 


600 




24 3 .48 


0.1550 


1730 


1.0 


95.6 


8.04 


5.5 


10 


-5 


2.2 


200 


600 


c s 


24 3 .48 


0.1570 


2110 


1.0 


92.4 


7.2 


5.5 


10" 


-6 


2.2 


400 


1200 



respectively. The resulting gluinoball propagator consists of connected and disconnected contribu- 
tions: 

disconnected connected 



The connected term can be used to extract the adjoint pion mass M a - n (the last term in Eq. ( |4.l| ) 
vanishes for Y = 75). OZI arguments [Q] suggest M^_ K oc nig for light gluinos. 

The disconnected propagators were computed by using the Stochastic Estimators Technique 
(SET) in the spin dilution variant to reduce the large variance. As it is the case in QCD, the 
disconnected diagrams are intrinsically noisier than the connected ones and dominate the level 
of noise in the total correlator. We performed tests on few configurations with up to Af = 40 noisy 
estimators in order to study the limit at which only the gauge noise dominates the statistical error of 
the disconnected correlator. The optimal number of estimators was finally fixed between 16 and 22 
estimates for all runs. In the case of a-r\' reasonable signal-to-noise ratio is obtained allowing the 
extraction of the mass from the mass fit. This was not possible for the a- fa; the latter has a nonzero 
overlap with the vacuum, and hence its correlator has a non-vanishing vacuum contribution (the 
last term in Eq. (p~T|)). This results in a much worse signal-to-noise ratio, and the effective mass 
could only be extracted at very short time separations. In the future we shall consider new variance 
reduction techniques for a more precise computation of the disconnected diagrams. 

ii) Gluino-glueballs. The gluino-glueballs (g-g) are spin-j colour singlet states of a gluon and 
a gluino. They are supposed to complete the Wess-Zumino supermultiplet of the adjoint mesons 
[Q]. The full correlator is built up from plaquettes connected by a gluino propagator line: 

c ">o = -zE I {<%*^M<m^w^»w<^<£ p ) ■ (4 - 2) 

x,y i,j,k,l 
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Gluino-Glue and r\' mass 
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Figure 1: The masses of the low lying bound states of the jV 
multiplied by ro as a function the squared adjoint pion mass. 



1 SU(2) SYM theory. Shown are the masses 



The above correlator is a matrix with two independent components in Dirac space: 

cf g (At)=C l (At)S a P+C n (At)yf . (4.3) 

We see agreement in the masses extracted from each component, and we choose the time antisym- 
metric component C\ to fit the masses. We apply APE smearing for the links and Jacobi smearing 
for the fermion fields in order to optimise the signal-to-noise ratio and to obtain an earlier plateau 
in the effective mass. 

iii) Glueballs. According to []7jj, ± glueballs are expected to be members of a second Wess- 
Zumino supermultiplet. Their study in SYM presents difficulties which closely resemble those 
encountered in glueball spectroscopy in QCD and, fortunately, can be overcome with the same 
type of techniques, namely, APE smearing with the variational method. Also in this case we use 
the simplest interpolating operator for the scalar glueball ++ built from single space-like plaquette. 

It turns out, however, that the present statistics is not enough to obtain a reliable determination 
of the glueball masses. Therefore, we are planning to increase the statistics. 

Results. The masses of the a-T]' and the gluino-glueball are displayed in Fig. [I] as a func- 
tion of the squared adjoint pion mass in units of the Sommer scale parameter. Both bound state 
masses appear to be characterized by a linear dependence on (roM K ) 2 . The gluino-glueball turns 
out to be appreciably heavier (50%) than the a-r\' . Runs with and without Stout-smearing give 
consistent results for the a-r]' mass, while a discrepancy is observed for the gluino-glueball, which 
can be interpreted as an 0(a) discretisation effect. The comparison of the two runs at (roM^) 2 ~ 
4 in (2fm) 3 and (3fm) 3 volumes reveals small finite volume effects. The linear extrapolation 
of the a-r\' mass to massless adjoint pion (including the four lightest points) gives: roM a .„i = 
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1.247(48) [499(20)MeV]. A rough estimate for the gluino-glueball gives r M g . g = 2.5 -3 [1000- 
1200 MeV]. 

5. Summary and conclusions 

New results on the spectrum of bound states have been obtained by means of the TS-PHMC 
algorithm with noisy correction. The new algorithm, with an improved gauge action and Stout- 
smeared links, allows to obtain a significantly better performance compared to the previously used 
TSMB algorithm. This allows us to simulate in significantly larger volumes of (2fm) 3 and (3fm) 3 . 
Our first results for the masses are within errors equal in these two cases implying that these vol- 
umes are large enough for the study of the particle spectrum. We shall systematically investigate 
the finite volume effects in future publications. 

According to the low-energy effective theory of [||] the a-rj' and the gluino-glueball belong to 
the same supermultiplet and therefore should be degenerate in the SUSY limit. However our pre- 
liminary results show a gluino-glueball mass systematically heavier than the a-r\' until the lightest 
simulated gluino mass in the weakly broken SUSY region. Whether this outcome is a discretisation 
artifact or a physical effect will become clear in future studies at finer lattice spacings. If the latter 
case applies, the interpolating gluino-glueball operator could have dominant overlap with a mem- 
ber of a higher supermultiplet. The complete identification of possible supermultiplets requires the 
inclusion of the masses of the glueballs and the scalar a-fo bound states. 

The computations were carried out on Blue Gene L/P and JuMP systems at JSC Julich, 
Opteron PC-cluster at RWTH Aachen and the ZIV PC-cluster of the university of Munster (Ger- 
many). 
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